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Abstract. - A system exhibiting multiple simultaneously broken symmetries offers the oppor- 
tunity to influence physical phenomena such as tunneling currents by means of external control 
parameters. Time-reversal symmetry and inversion symmetry are both absent in ferromag- 
netic metals with substantial spin-orbit coupling. We here study transport of spin in a system 
consisting of two ferromagnets with spin-orbit coupling separated by an insulating tunneling 
junction. A persistent spin-current across the junction is found, which can be controlled in a 
well-defined manner by external magnetic and electric fields. The behavior of the spin-current 
for important geometries and limits is studied. 



o 



> 

o 



O 



l/"") \ Due to the increasing interest in the field of spintronics in recent years [1], the idea of 

utilizing the spin degree of freedom in electronic devices has triggered an extensive response 
in many scientific communities. The spin-Hall effect is arguably the research area which has 
received most focus in this context, with substantial effort being put into theoretical con- 
siderations [2] as well as experimental observations [3] . As the ambition of spintronics is to 

C^ ■ make use of the spin degree of freedom rather than electrical charge, investigations of mech- 
anisms that offer ways of controlling spin-currents are of great interest. Ferromagnetism and 
spin-orbit coupling are physical properties of a system that crucially influence the behavior of 
spins present in that system. For instance, the presence of spin-orbit coupling is highly im- 

S ' portant when considering ferromagnetic semiconductors [4]. Such heterostructures have been 
proposed as devices for obtaining controllable spin injection and manipulating single electron 
spins by means of external electrical fields, making them a central topic of semiconductor 
spintronics [5]. In ferromagnetic metals, spin-orbit coupling is ordinarily significantly smaller 
than for semiconductors due to the bandstructure. However, the presence of a spin-orbit 
coupling in ferromagnets could lead to new effects in terms of quantum transport. 

Studies of tunneling between ferromagnets have uncovered interesting physical effects 
[7,15]. Nogueira et al. predicted [15] that a dissipationless spin-current should be estab- 
lished across the junction of two Heisenberg ferromagnets, and that the spin-current was 
maximal in the special case of tunneling between planar ferromagnets. Also, there has been 
investigations of what kind of impact spin-orbit coupling constitutes on tunneling currents 
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in various contexts, e.g. for noncentrosymmetric superconductors [8], and two-dimensional 
electron gases coupled to ferromagnets [9] . Broken time reversal- and inversion-symmetry are 
interesting properties of a system with regard to quantum transport of spin and charge, and 
the exploitment of such asymmetries has given rise to several devices in recent years. For 
instance, the broken 0(3) symmetry exhibited by ferromagnets has a broad range of pos- 
sible applications. This has led to spin current induced magnetization switching [10], and 
suggestions have been made for more exotic devices such as spin-torque transistors [11] and 
spin-batteries [12]. It has also led to investigations into such phenomena as spin-Hall effect in 
paramagnetic metals [14], spin-pumping from ferromagnets into metals, enhanced damping of 
spins when spins are pumped from one ferromagnet to another through a metallic sample [13], 
and the mentioned spin Josephson effects in ferromagnet/ferromagnet tunneling junctions [15]. 

In this Letter, we investigate the spin-current that arises over a tunneling junction sep- 
arating two ferromagnetic metals with substantial spin-orbit coupling. It is found that the 
total current consists of three terms; one due to a twist in magnetization across the junction 
(in agreement with the result of ref. [15]), one term originating from the spin-orbit inter- 
actions in the system, and finally an interesting mixed term that stems from an interplay 
between the ferromagnetism and spin-orbit coupling. After deriving the expression for the 
spin-current between Heisenberg ferromagnets with substantial spin-orbit coupling, we con- 
sider important tunneling geometries and physical limits of our generally valid results. Finally, 
we make suggestions concerning the detection of the predicted spin-current. Our results indi- 
cate how spin transport between systems exhibiting both magnetism and spin-orbit coupling 
can be controlled by external fields, and should therefore be of considerable interest in terms 
of spintronics. 

Our system consists of two Heisenberg ferromagnets with substantial spin-orbit coupling, 
separated by a thin insulating barrier which is assumed to be spin-inactive. This is shown in 
fig. m A proper tunneling Hamiltonian for this purpose is Ht = X)kpo-(-^' k P c k(T C W "+" h.c.), 
where {c^Cko-} and {d^ , dko-} are creation and annihilation operators for an electron with 
momentum k and spin a on the right and left side of the junction, respectively, while T^ p is 
the spin-independent tunneling matrix element. In k-space, the Hamiltonian describing the 
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Fig. 1 - Our model consisting of two ferromagnetic metals with spin-orbit coupling separated by 
a thin insulating barrier. The magnetization m and electrical field E are allowed to point in any 
direction so that our results are generally valid, while special cases such as planar magnetization etc. 
are easily obtained by applying the proper limits to the general expressions. 
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in which e^ is the kinetic energy of the electrons, J is the ferromagnetic coupling constant, 
N is the number of particles in the system, 77 (k) is a geometrical structure factor, while Sk = 
(1/2) Yl a c iia a 'apCkP 1S t ne s P m operator. As we later adopt the mean-field approximation, 
m = (m x , Tn y , m z ) will denote the magnetization of the system. 

The spin-orbit interactions are accounted for by a Rashba Hamiltonian 

#S-0 = -$>k[£( Wxk )-^k, (2) 

k 

where </j k = [ckj,c k J T , E = — W is the electrical field felt by the electrons and a = 
(cii Gii a 3) m which Oi are Pauli matrices, while the parameter £ is material-dependent. From 
now on, the notation £(E x k) = Bk = (-B x ,k> B y k, B z k) will be used. In general, the elec- 
tromagnetic potential V consists of two parts Vnt and Vo Xt (see e.g. ref. [5] for a detailed 
discussion of the spin-orbit Hamiltonian) . The crystal potential of the material is represented 
by Vint, and only gives rise to a spin-orbit coupling if inversion symmetry is broken in the 
crystal structure. Asymmetries such as impurities and local confinements of electrons are 
included in V cx t, as well as any external electrical field. Note that any lack of crystal inver- 
sion symmetry results in a so-called Dresselhaus term in the Hamiltonian, which is present 
in the absence of any impurities and confinement potentials. In the following, we focus on 
the spin-orbit coupling resulting from Vext , thus considering any symmetry-breaking electrical 
field that arises from charged impurities or which is applied externally. Assuming that the 
crystal structure respects inversion symmetry, a Dresselhaus term [6] is nevertheless easily 
included in the Hamiltonian by performing the substitution (E x k) • <x — > [(E x k) + P(k)] • <r, 
where £>(k) = -D(-k). 

We now proceed to calculate the spin-current that is generated across the junction as a 
result of tunneling. Note that in our model, the magnetization vector and electrical field are 
allowed to point in arbitrary directions. In this way, the obtained result for the spin-current 
will be generally valid and special cases, e.g. thin films, are easily obtained by taking the 
appropriate limits in the final result. 

In the mean-field approximation, the Hamiltonian for the right side of the junction can be 
written as H = -ffpM + ^s-o 1 which in a compact form yields 

where £k<r = £k — &{Cz,r — B z .y) and Ho is an irrelevant constant. The FM order parameters are 
£r, = 2 Jrj(0)(m x ^ — iniy^B.) and Cz,R. = 2 Jrj(0)m z ^ and /Ck = B x ^ — i-B y ,k- For convenience, 
we from now on write £ = ICIe 1 ^ and /Ck = |/Ck|e l6lk . The Hamiltonian for the left side of the 
junction is obtained from eq. (0) simply by the doing the replacements k — > p and R — > L. 

In order to obtain the expressions for the spin- and charge- tunneling currents, it is neces- 
sary to calculate the Green functions. These are given by the matrix Gk{^ n ) — (~ i w n+-4k) , 
where Ak is the matrix in eq. (0). Explicitly, we have that 

Above, ui n = 2(n + l)ir/f3, n = 0, 1, 2 ... is the fermionic Matsubara frequency and (3 denotes 
inverse temperature. Introducing X k (io; n ) — (e k -f ~ iw„)(£ k | — iw n ) — |Cr — /C k | 2 , the normal 
and anomalous Green's functions are 

G k CT (i^„) = (e k ,- CT - kj n )/X k (iu> n ), F k iT (i^„) = F^(iw n ) = (Cr - /C k )/X k (iu;„). (5) 
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Defining a spin-current is not as straight-forward as defining a charge-current [16]. Specifi- 
cally, the conventional definition of a spin-current given as spin multiplied with velocity suffers 
from severe flaws in systems where spin is not a conserved quantity. In this Letter, we de- 
fine the spin-current across the junction as I s (i) = (dS(i)/di) where dS/di = i[i7x,S]. It is 
then clear that the concept of a spin-current in this context refers to the rate at which the 
spin- vector S on one side of the junction changes as a result of tunneling across the junction. 
We choose the right side, denoting the physical parameters with labels k and R. In this way, 
we avoid non-physical interpretations of the spin-current in terms of real spin transport as 
we only calculate the contribution to dS/di from the tunneling Hamiltonian instead of the 
entire Hamiltonian H. Should we have chosen the latter approach, one would run the risk of 
obtaining a non-zero spin-current due to e.g. local spin-flip processes which are obviously not 
relevant in terms of real spin transport across the junction. 

The expression for I (t) is established by first considering the generalized number op- 
erator N a p — X)k c ka c k/3- This operator changes with time due to tunneling according 
to N a p — i[H T ,N a p], which in the interaction picture representation becomes N a p(t) = 
— iSkp(^kp c ka^p/3 el * ey — h.c). The voltage drop across the junction is given by the differ- 
ence in chemical potential on each side, i.e. eV = /xr — /il- The Matsubara formalism then 
dictates that the spin-current across the junction is 

!*(*) = 5 £*<*<#*(')>. ( 6 ) 

a/3 

where the expectation value of the time derivative of the transport operator is calculated in the 
linear response regime by means of the Kubo formula (N a p(t)) = —' 1 J_ 00 dt'([N a p(t), ifr (£')])• 
Consider now the ^-component of the spin-current in particular, which can be written 
as if = 3m{$(— eV)}. The Matsubara function $(— eV) is found by performing analytical 
continuation ip„ — * —eV + i0 + on $(ip„), where 

*(iPn)=i Yl ^^(GrK„)G7(i Wm -ip„) + F k -^ CT (i^„)^~ CT (i^ m -iPn))- (7) 

Here, p n = 2nir/f3, n = 0, 1, 2 ... is the bosonic Matsubara frequency. Inserting the Green's 
functions from eq. (|g) into eq. ([?]), one finds that a persistent spin-current is established 
across the tunneling junction. For zero applied voltage, we obtain 



s = y^ \ T *p\ J k P t |CrCl i sin A0 + |/Ck/C | sin M 



k P 



^k P = y^ afi 



a=± 
p=± 



- |£k<L|sin(0 k -0 L ) - |£ P CR|sin(<?!>R.-0p) , (8a) 

n(e k + aryk) - n(e p + f3j p ) 
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(8b) 



In eqs. (§), A0 kp = k - p , A^ = 0R - <i>U while 7^ = (Q tR - B zM ) 2 + \( R - /C k | 2 and 
n(e) denotes the Fermi distribution. In the above expressions, we have implicitly associated 
the right side R with the momentum label k and L with p for more concise notation, such 
that e.g. i? 2jk = Bf k . The spin-current described in eq. (||) can be controlled by adjusting 
the relative orientation of the magnetization vectors on each side of the junction, i.e. A</>, 
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and also responds to a change in direction of the applied electric fields. The presence of an 
external magnetic field H^ would control the orientation of the internal magnetization m;. 
Alternatively, one may also use exchange biasing to an anti-ferromagnet in order to lock the 
magnetization direction. Consequently, the spin-current can be manipulated by the external 
control parameters {Hi,Ei} in a well-defined manner. This observation is highly suggestive 
in terms of novel nanotechnological devices. 

We stress that eq. m) is non-zero in the general case, since 7k ^ — 7— k and #_k — 9^ + it. 
Moreover, eq. (H) is valid for any orientation of both m and E on each side of the junction, 
and a number of interesting special cases can now easily be considered simply by applying the 
appropriate limits to this general expression. 

Consider first the limit where ferromagnetism is absent, such that the tunneling occurs 
between two bulk materials with spin-orbit coupling. Applying m — > to eq. (@), it is readily 
seen that the spin-current vanishes for any orientation of the electrical fields. Intuitively, one 
can understand this by considering the band structure of the electrons and the corresponding 
density of states N(e) when only spin-orbit coupling is present, as shown in fig. 0. Since the 
density of states is equal for f and j spins, one type of spin is not preferred compared to the 
other with regard to tunneling, resulting in a net spin-current of zero. Formally, the vanishing 
of the spin-current can be understood by replacing the momentum summation with integration 
over energy, i.e. X)k P ~^ / I ^R^hN^e^N^ei,). When m — > 0, eq. (|j) dictates that 



/3=± 
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deR, Q de L ,/3A r R(eR,a)A r L(eL, / 3) 
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(9) 



Since the density of states for the f- and ^-populations are equal in the individual subsystems, 
i.e. N^(e) — N^-(e) = N(e), the integrand of eq. (g) becomes spin-independent such that 
the summation over a and (3 yields zero. Thus, no spin-current will exist at eV — over a 
tunneling barrier separating two systems with spin-orbit coupling alone. In the general case 
where both ferromagnetism and spin-orbit coupling are present, the density of states at, say, 
Fermi level are different, leading to a persistent spin-current across the junction due to the 
difference between N'(e) and N^(e). 



a) 





Fig. 2 - Schematic illustration of the energy-bands for a) a system with spin-orbit coupling, b) a 
system with ferromagnetic ordering, and c) a system exhibiting both of the aforementioned properties. 
Since the density of states N a (£fc CT ) is proportional to (deka/dk) -1 , we see that a difference between 
N (sfe CT ) and iV (sfc CT ) is zero in a), while the density of states differ for the |- and .[-populations in b) 
and c). Thus, a persistent spin-current will only occur for tunneling between systems corresponding 
to b) and c). 
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We now consider a special case where the bulk structures indicated in fig. [j] are reduced to 
two thin-film ferromagnets in the presence of electrical fields that are perpendicular to each 
other, say El = (-El, 0, 0) and Er = (0, E-r, 0), as shown in fig. Em) and b). In this case, we 
have chosen an in-plane magnetization for each of the thin-films. Solving specifically for fig. 
Ha), it is seen that niL = (OjTO^l, jti z ,l) and rriR = (m x ,R, 0, TO z ,r) . Furthermore, assume 
that the electrons are restricted from moving in the "thin" dimension, i.e. p = (0,p y ,p z ) and 
k = (k x , 0,k z ). In this case, eq. (pf) reduces to the form 



■Jf = 4sgn(TO yjL ) + Y^ ■Ji,kpSgn(p*) 
k P 



(10) 



where ii,kp ^ Ii,-k,-p', details of the calculations for the setup in fig. 0a) and b) will be given 
in a more comprehensive paper [18]. 

From these observations, we can draw the following conclusions: whereas the spin-current 
is zero for the system in Fig. 0a) if only spin-orbit coupling is considered, it is non-zero 
when only ferromagnetism is taken into account. However, in the general case where both 
ferromagnetism and spin-orbit coupling are included, an additional term in the spin-current is 
induced compared to the pure ferromagnetic case. Accordingly, there is an interplay between 
the magnetic order and the Rashba-interaction that produces a spin-current which is more 
than just the sum of the individual contributions. 
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Fig. 3 - Tunneling between planar ferromagnets in the presence of externally applied electrical fields 
El and Er that destroy inversion symmetry and induce a spin-orbit coupling. 



Detection of the induced spin-currents would be challenging, although recent studies sug- 
gest feasible methods of measuring such quantities. For instance, the authors of ref. [20] pro- 
pose a spin-mechanical device which exploits nanomechanical torque for detection and control 
of a spin-current. Similarly, a setup coupling the electron spin to the mechanical motion of 
a nanomechanical system is proposed in [19]. The latter method employs the strain-induced 
spin-orbit interaction of electrons in a narrow gap semiconductor. In rcf. [23], it was demon- 
strated that a steady-state magnetic-moment current, i.e. spin-current, will induce a static 
electric field. This fact may be suggestive in terms of detection [21,22], and could be useful 
to observe the novel effects predicted in this Letter. 

In summary, we have derived an expression for a dissipationless spin-current that arises 
in the junction between two Heisenberg ferromagnets with spin-orbit coupling. We have 
shown that the spin-current is driven by terms originating from both the ferromagnetic phase 
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difference, in agreement with the result of ref. [15], and the presence of spin-orbit coupling 
itself. In addition, it was found that the simultaneous breaking of time-reversal and inversion 
symmetry fosters an interplay between fcrromagnetism and spin-orbit coupling in the spin- 
current. Availing oneself of external magnetic and electric fields, our expressions show that 
the spin-current can be tuned in a well-defined manner. These results are of significance in the 
field of spintronics in terms of quantum transport, and offer insight into how the spin-current 
behaves for nanostructures exhibiting both ferromagnctism and spin-orbit coupling. 
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